Kinetics of a frictional granular motor 
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Within the framework of a Boltzmann-Lorentz equation, we analyze the dynamics of a granular rotor 
immersed in a bath of thermalized particles in the presence of a frictional torque on the axis. In numerical 
simulations of the equation, we observe two scaling regimes at low and high bath temperatures. In the 
large friction limit, we obtain the exact solution of a model corresponding to asymptotic behavior of the 
Boltzmann-Lorentz equation. In the limit of large rotor mass and small friction, we derive a Fokker-Planck 
equation for which the exact solution is also obtained. 
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Recently, Eshuis et al. 0], inspired by Smoluchowski's 
gedankenexperiment, constructed a macroscopic rotational 
motor consisting of four vanes immersed in a granular gas. 
When a soft coating was applied to one side of each vane, 
a ratchet effect was observed above a critical granular tem- 
perature. While this was the first experimental realiza- 
tion of a granular motor, similar Brownian ratchets exist in 
many diverse applications, e.g., photovoltaic devices and 
biological motors; See (2). All of these motors share the 
common features of non-equilibrium conditions and spa- 
tial symmetry breaking. Several recent theoretical stud- 
ies of idealized models of granular motors, which use a 
Boltzmann-Lorentz description J3j-[H1» show that the motor 
effect is particularly pronounced when the device is con- 
structed from two different materials, as was the case in 
the recent experiment Q. The existing theories, however, 
predict a motor effect for any temperature of the granular 
gas while in the experiment the phenomenon is only ob- 
served if the bath temperature is sufficiently large. The 
presence of friction is at the origin of this difference and 
it is therefore highly desirable to incorporate it in the theo- 
retical schemes. 

Several studies have considered the effect of friction on 
Brownian motion. The pioneering work of de Gennes |9) 
was motivated by the motion of a coin on a horizontally 
vibrating plate[10| and a liquid droplet on non-wettable 
surfaces subjected to an asymmetric lateral vibration; he 
obtained the stationary velocity distribution, as well as the 
velocity correlation function. A notable result, obtained 
independently by Hayakawa lUTl . is that the Einstein rela- 
tion no longer holds in the presence of dry friction. More 
recently, Touchette and coworkers lfl"2i4l"4ll obtained a so- 
lution of a model with dry friction and viscous damping. 

In this Letter, we investigate the kinetic properties of 
a heterogeneous granular rotor in the presence of dry 
(Coulomb) friction and we provide exact solutions in the 
limit of large motor mass (Brownian limit) and also in the 
limit of strong friction for arbitrary mass . 

We consider a two-dimensional system where a hetero- 



geneous chiral rotor with moment of inertia /, mass M and 
the length L immersed in a granular gas at density p with a 
velocity distribution (f)(v) characterized by a granular tem- 
perature T. The motor is composed of two materials with 
coefficients of restitution a + and a_ . Collisions of the bath 
particles, of mass m, with the former (latter) exert a posi- 
tive (negative) torque. (The experimental setup employed 
by Eshuis et al., where one has a translational invariance in 
the third dimension, can also be described by the model). 




FIG. 1 : The chiral rotor is constructed of materials with coeffi- 
cients of restitution a+ (red) and a_ (blue). When immersed in 
a bath of thermalized particles, it acquires a net rotation around 
its axis. The model can be extended to include additional vanes 
(shown in gray) as used in the experiment of Eshuis et al. (T). 

To understand the qualitative effect of the friction, it 
is useful to consider the torques acting on the rotor: The 
motor torque is the average of all changes of angular mo- 
mentum resulting from collisions between the rotor and the 
bath particles. In the large rotor mass ratio (M/m >> 1), 
one finds that the motor torque is (a + — a_)^pL 2 T lfT5l . 
which vanishes for a symmetric rotor (a + = a_) and the 
dynamic friction is assumed constant and opposed to the 
motion -T s a(n) (a(±x) = ±1, \x\ ^ 0; <r(0) = 0). 
(We neglect static friction in this study, but it will be con- 
sidered in a future publication). We can define a crossover 
temperature, T c = pL 2^' + s _ a y that separates a regime 
dominated by the motor effect (T > T c ) and one domi- 
nated by friction (T < T c ). This mechanical point of view 
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defining the two kinetic regimes must be expanded by ac- 
counting for the fact that the motor torque is a fluctuating 
quantity. 

Assuming the bath is homogeneous, and in the presence 
of dry friction, the Boltzmann-Lorentz (BL) equation of 
the angular velocity distribution (AVD), /(O), can be ex- 
pressed as 



at 



/(«;*) 



I dti 



f(n;t) = J[f,cf>] (1) 



where the frictional torque T s and J[f, (ft] is the BL colli- 
sion operator |[T5l 
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with O**, -y** denoting the pre-collisional velocities, 9{x) 
is the Heaviside function and pv(fl) is the total colli- 
sion rate for a motor rotating at the angular velocity Q, 



pv(Q) = pJZo dv I l/2 dX \ v ~ Afi|0(t;). (For a 3D- 
motor, Eq{T]is unchanged if p = nLi, n being the 3D 
number density and L\ the perpendicular linear size of the 
vane.) For convenience we introduce the dimensionless 
variables Q* = LJ^Sl, T* = and v* = 



T uu ' *■ « pL 2 T " " V T ' 

With this choice of units, the steady state angular velocity 
depends only on a±, M/m and T*. Also, since the BL de- 
scription neglects recollisions with the bath particles and is 
invalid when M « m, we will assume in the following 
that M > m. 

Simulations. A complete analytical treatment of the 
Boltzmann equation is a very difficult task, even in the 
steady state. We have performed numerical simulations of 
the BL equation. (See |fl6l for details. We note that both 
in simulation and experiment the sampling time must be 
smaller than the mean stopping time to avoid biasing the 
AVD). Fig. [2] shows the dimensionless mean angular ve- 
locity O* of an asymmetric rotor (a+ = 1, a_ = 0) as 
a function of a dimensionless dry friction T* for different 
mass ratios M/m = 1,2, 10, 20. As expected, the mean 
angular velocity decreases as the friction increases. More- 
over, two scaling regimes are observed: 
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In experiments, the most easily varied parameter is the 
bath temperature T, and in terms of the original variables, 
the two scaling regimes correspond to (Q) oc T 3 / 2 and 
(Q) ~ T 1 / 2 at low and high temperatures, respectively. 
(The inset of Fig.|2]shows a log-log plot of the mean angu- 
lar velocity (O) as a function of bath temperature T.). 

Limiting Behaviors. We now show that in appropriate 
limits, the physics can be described by simpler models for 




FIG. 2: (Color online) Dimensionless mean angular velocity 
(ft*) of an asymmetric rotor versus 1^: a+ = 1, a— = Ofor 
different mass ratios M/m — 1, 2, 10, 20, top to bottom obtained 
from numerical simulations of the BL equation. The dashed 
and dotted curves correspond to the solutions of the Independent 
Kick model and the Fokker-Planck equation, respectively. The 
inset shows (fi) as a function of the bath granular temperature, 
T, for T s = 0.2 and M/m = 5,10,20. At low temperature 
< Q >~ T 3/2 and at high temperature < fi >~ T 1/2 . The 
dashed lines are visual guidance for this scaling behavior. 



which analytical solutions can be obtained. The analy- 
sis is facilitated by considering two time-scales: the av- 
erage inter-collision time between bath particles and the 
motor, r c 



r ^—F (for a Gaussian bath dis- 
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tribution) and the mean stopping time in the presence of 
friction, t s = |P where ft is the average angular ve- 
locity after a collision. When the dimensionless friction 
r* is small, Q* ~ (a + — a_), while for V* » 1, 
Q* ~ 1I ^(a + — a_), which gives 



t s a+ - 



1, 
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Whenever slip-stick behavior is present the AVD con- 
tains a regular part and a delta singularity at = corre- 
sponding to a stationary rotor before the next collision with 
a bath particle: 

/*(fi*)=7/fl(n*) + (l-7W) (5) 

where J dO,* /^(O*) = 1 and 7 is a constant that can be 
determined from conservation of the probability current at 
n* = OOl: 
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with C 



2/^(0) / J_ dy\y\<ft(y) is a numerical constant 
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When r c > > t s the dry friction stops the motor before 
the next collision, and the motor essentially evolves by fol- 
lowing a sequence of stick- slip motions. Most of the time, 
the rotor is at rest and the singular contribution is dominant, 
7 ~ 1/r*. This regime can be described by the Indepen- 
dent Kick model introduced below. 

Conversely, when t c << r s , collisions are so frequent 
that sliding dominates the motor dynamics. For all prac- 
tical purposes, the rotor never stops and (1 — 7) ~ T*. 
In this case, the dynamics is well described by a Fokker- 
Planck equation for M/m >> 1. 

Independent Kick Model. The mean angular velocity is 
the average over all collisions: 
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angular velocity after a collision, which is given by 



and ft is the 
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Integrating over time and x, and by using a Gaussian bath 
distribution (analytical expressions can be obtained for ar- 
bitrary bath distributions, but they are not presented here), 
the dimensionless mean angular velocity is given by 
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where £ = 12 j]~- Comparison of the exact expression of 
the Independent Kick model with the simulation results of 
the BL equation in Fig. [2] shows that it has the correct di- 
mensionless dependence 1 /T* . 

Finally, the regular part of the AVD can be obtained from 
the inverse Fourier transform of the characteristic func- 
tion {e lka ). For a Gaussian bath distribution, by using the 
change of variable z = 2x/L, one obtains 



1 



9(fl*) / dze 8 « 2 d+°+> 2 



+0(-ft*) / dze 



(i+;»-)-s)'- 



(10) 



Therefore, the Independent Kick model provides non- 
trivial angular velocity distributions which are the sum of 
a regular function with a non-Gaussian shape, even in the 
Brownian limit, with an amplitude decreasing as 1 /Y* and 
a singular contribution. Fig. [5] a displays regular angular 
velocity distributions of the BL model for different fric- 
tional torques. As expected the amplitude decreases as the 
friction increases. Fig. [3]b shows the rescaled distributions 
T*7/^(0*) versus ft*. One observes that for Y* > 2, 
the curves converge to the AVD of the Independent Kick 



model, Eq. ( 10 1 confirming that it describes very accurately 
all kinetic properties of the BL model for moderate to large 
friction. 




FIG. 3: (Color online) (a) Regular part of the AVD, 7/£(fi*), for 
various frictional torques T* = 1.22, 2.72, 6.05, 13.6, 30.0 (From 
top to bottom) obtained from numerical simulations of the BL 
equation. M/m = 10 with other parameters the same as in Fig. 
[2] (b) rescaled T* sl f^(Q*) for the same values of the friction 
torque (from bottom to top). The upper curve corresponds to 
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Eq[l0] (c) AVD (full curves) and FP solutions (Eq_ 

dashed curves) for M/m = 20 and T* = 0.01, 0.033, 0.074, 0.1 
(d) Log-linear plot of for the same values of the friction 

torque. 



The amplitude of the singular contribution of 
is given by 1 — 7 = 1 — (e zkn *) k =o = 1 _ 
^ 2+ 8g"r* a ~^ m (1 + £)• Consequently, the model remains 
valid if 7 < 1, namely if the frictional torque is sufficiently 
large. For moderate friction, the Independent Kick model 
overestimates the mean angular velocity given by the BL 
equation (see Fig.|2]l. 

Brownian Limit and the Fokker-Planck Equation. We 
now consider the opposite limit when the stopping time 
t s is much larger than the intercollision time r c and when 
M » m. For homogeneous granular motors where the 
mean angular velocity goes to zero in the Brownian limit 
l3l[T2l, the Fokker-Planck can be derived by performing a 
Kramers-Moyal (KM) expansion of the BL integral oper- 
ator. For heterogeneous granular motors, where the mean 
angular velocity remains finite, we introduced a KM-like 
method providing the exact Brownian limit[15|, but finite 
mass corrections cannot be easily included. We propose 
here an approach that consists of performing a complete 
series expansion in terms of derivatives of /(f)) and a com- 
plete resummation of the coefficients allowing the BL op- 
erator to be expressed as 
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where the generating function 
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g(tt,a) is given by 



g(Q, a) = 2p 
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Truncating the BL operator, Eq. (ill, at second-order and 



adding the frictional torque leads to the following Fokker- 
Planck equation 
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in which all finite -mass corrections are incorporated, and 
where deviations from a Gaussian distribution are present 
for large finite masses. Recalling that the g n (Q) are pro- 
portional to p, we see that for a given frictional torque 
r s , increasing the bath density reduces the effect of fric- 
tion. The corresponding Langevin equation features a mo- 
tor torque with a non-linear dependence on f2 and a colored 
noise. 

The steady state solution of Eq. ( fT3) l is 
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du 



gi(tt) + T s a{u) 
92(u) 
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where Cr 3 is obtained from the normalization condition 
J dQ,f(Q) = 1. This result clearly shows that, even in 
the absence of friction, the AVD is non-Gaussian for finite 
mass ratios. Moreover, the signature of the presence of 
friction is a cusp of /(f)) at O, = 0. This feature, also ob- 
tained in the Independent Kick model, is always observed 
in simulation results of the BL equation. 

In the Brownian limit <?i(J2) ~ g[(Cl)(Q — (l) and 
g2(yi) = 2T g g[(Q), where T g (the rotor granular temper- 
ature which is lower than the bath temperature, T) and Q 
are given by the Kramers-Moyal expansion (Eqs. 10 and 
12 in Ref. 03)). This finally gives a stationary distribution, 
at the lowest order in m/M, 



/(fi) = Cexp 



i(n -nf ir s \n\ 



2T n 



dm 
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where C is the normalization constant. Whereas one ob- 
serves a Gaussian decay of the AVD at large velocity, /(fi) 
decreases exponentially for small and intermediate angular 
velocities, due to the presence of friction. The analytical 
expression (EqfT5j) provides an accurate description of the 
AVD obtained by simulation (see Fig(3]c and d). A first- 
order expansion of the dimensionless mean angular veloc- 
ity can be easily obtained 



where a x and a 2 are positive constants. The dominant 
term corresponds to the first scaling law obtained in sim- 
ulation results of the BL equation. More generally, the 
prediction of the Fokker-Planck equation (dotted curve) is 
compared to simulation results of the BL equation as a 
function of the frictional torque T* for M/m = 10 (see 
Fig.|2|. As expected, the Fokker-Planck description is ac- 
curate at low and moderate friction, whereas the Indepen- 
dent Kick model is asymptotically exact in the large frac- 
tion limit. Even in the small region of moderate frictions 
r*, the two models slightly overestimate the mean angular 
velocity. We note that Touchette et al. obtained a time- 
dependent solution of the FP equation of a particle subject 
to solid friction and viscous damping lfT2l . corresponding 
to a + = a_ in the present model. This approach can be 
extended to obtain full solution of the granular Brownian 
motor. 

In their experimental study [1], Eshuis et al. observed 
that the mean angular velocity increases as (S — S C ) 1A 
where S is the dimensionless shaking strength and S c is 
a threshold value. The relation between S — S c and T is 
not clear, but our analysis suggests two scaling laws could 
be observed in experiments corresponding to the small and 
large friction limits. Further experiments should clarify 
this. 

We thank Olivier Benichou for helpful discussions. 
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